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Abstract. We determine the weak limit of the distribution of the random variables "height" and "range" on the 
set of p-watermelons without wall restriction as the number of steps tends to infinity. Additionally, we provide 
asymptotics for the moments of the random variable "height" . 



1. Introduction 

The model of vicious walkers was originally introduced by Fisher [6] as a model for wetting and melting processes. 
In general, the vicious walkers model is concerned with p random walkers on a <i-dimensional lattice. In the lock 
step model, at each time step all of the walkers move one step in one of the allowed directions, such that at no 
time any two random walkers share the same lattice point. 

A configuration that attracted much interest amongst mathematical physicists and combinatorialists is the wa- 
termelon configuration, which is the model underlying this paper (see Figure [T] for an example). This configuration 
can be studied with or without the presence of an impenetrable wall. By tracing the paths of the vicious walkers 
through the lattice we can identify the (probabilistic) vicious walkers model with certain sets of non-intersecting 
lattice paths. It is exactly this equivalent point of view that we adopt in this paper. We proceed with a precise 
definition. A p-watermelon of length 2n is a set of p lattice paths in Z 2 satisfying the following conditions: 

• the z-th path starts at position (0, 2i) and ends at (2n, 2i), 

• the paths consist of steps from the set {(1,1), (1,-1)} only and 

• the paths are non-intersecting, that is, at no time any two path share the same lattice point. 

An example of a 4-watermelon of length 16 is shown in Figure [T] (for the moment, the dashed lines and the labels 
should be ignored). 

Since its introduction, the vicious walkers model has been studied in numerous papers. While early results 
mostly analyse the vicious walkers model in the continuum limit, there are nowadays many results for certain 
configurations directly based on the lattice path description given above. With the increasing number of results it 
became clear that vicious walkers are very important objects in mathematical areas far beyond its original scope. 
For example, Guttmann, Owczarek and Viennot [10] related the star and watermelon configurations to the theory 
of Young tableaux and integer partitions. Later, Krattenthaler, Guttmann and Viennot [16] proved new, exact as 
well as asymptotic, results for the number of certain configurations of vicious walkers. 

The vicious walkers model is also very closely related to random matrix theory, as can be seen from articles by, 
e.g., Baik [1], Johansson [11] and Nagao and Forrester [18]. Later, Katori and Tanemura [13] and Gillet [9] studied 
the diffusion scaling limit of certain configurations of vicious walkers, namely stars and watermelons, respectively. 

In 2003, Bonichon and Mosbah [2] presented an algorithm for uniform random generation of watermelons, which 
relies on the counting results by Krattenthaler, Guttmann and Viennot [16]. Amongst other things, Bonichon and 
Mosbah studied the parameter height on the set of watermelons (with and without wall) . 

In this paper we rigorously analyse the following two parameters on the set of p-watermelons: 

• The height of a watermelon is the maximum ordinate reached by its top most branch. 

• The range of a watermelon is the difference of the maximum of its top most branch and the minimum of 
its bottom most branch (the depth of the watermelon) . 

The 4-watermelon depicted in Figure [T] has the height 11 and the range 11 + 4 = 15. 
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Figure 1. Example of a 4-watermelon of length 16 without wall, height 11, depth 4 and range 15 

Katori et. al. [12] and Schehr et. al. [19] studied the parameter "height" in the continuous limit, and recovered 
the leading terms for some of the asymptotics proved in [4,5]. Additionally, Schehr et. al. gave some arguments 
concerning the behaviour of the parameter "height" as the number of walkers tends to infinity. 

Now, consider the set of p-watermelons of length 2n, endowed with the uniform probability measure. We 
can then speak of the random variables "height", denoted by H ntP , and "range", denoted by R n , P , on this set. 
We determine the weak limits of H n p and R n , P as the number n of steps tends to infinity (see Theorem [T] and 
Theorem [31 respectively). Additionally, we determine asymptotics for the moments of H n ^ p . In particular, we 
prove that the s-th moment of the random variable "height" behaves like K s n s / 2 + r s n < ' s ~ 1 ^ 2 + (n s / 2_1 ) for some 
explicit numbers k s and r s , see Theorem [2] 

Techniques similar to those applied in this paper can also be used to analyse the random variable height on the 
set of p- watermelons under the presence of an impenetrable wall. For details we refer to [5] . 

The paper is organised as follows. The next section contains some well known results that are needed in the 
subsequent sections. In Section [3J we consider the random variable "height" , and we determine the weak limit 
as well as asymptotics for all moments. In the last section, we determine the weak limit of the random variable 
"range" . 

2. Preliminaries 

In this section we collect several results which will be needed in the two subsequent sections. All these results are 
either well known in the literature and/or can easily be derived by means of standard techniques. We, therefore, 
remain very brief, give only a few comments on the proofs and in each case refer to the corresponding literature 
for details. 

We start with an exact enumeration result for the total number of watermelons confined to a horizontal strip. 
Lemma 1. The number m^ h k of p-watermelons without wall, length 2n, height < h and depth > ~k is given by 

m (rt = det (VlY ^ \-( ^ 

n ' h ' k o<i,j<pyf-^\\n + £(h+k) + i- jj \n + i(h + k) + h - i - j 
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The total number m„ of p-watermelons is given by 



= det 

o<i,j<p \ \n + i — j 



2n 



This lemma follows immediately from the well-known Lindstrom-Gessel-Viennot formula (see [8, Corollary 3] 
or [17, Lemma 1]), together with an iterated reflection principle. 



Remark 1. Since any p- watermelon without wall and length 2n has depth > — n — 1, we see that the number of 
watermelons with height < h and no restriction on the depth is given by 

^vT/i n+i* For the sake of convenience, 
this quantity will also be denoted by m^ h . In this special case, the determinantal expression above simplifies to 



m«= det 2n )- 2U 

n ' h o<i,j< P \ \n + i — jj \n + h — i — j 



Lemma 2. We have 

•2\® /2nWCi 



m( n ] = 1-1 



n J \n 



i=0 / 



as n ^ oo. 

» 



Proof (Sketch). The result is established from the closed form expression for m„ , viz 

»- de, (( 2 » )W 2 1Ttp- <2 " n ' x2N 

u \ \ r> -I- ? — til \ n I \ 



m " "oim n |i-jjJ"i n J \JLL*- (2n)! V(n + 0! 
For details on the evaluation of this (and many more) determinant, we refer to [15]. □ 
Lemma 3. For \m — z\ < n 5 / 8 , z bounded, and arbitrary N > 1 we have the asymptotic expansion 

In \ 47V+1 

In I U\ 



m U+m-J _ -rn 2 /n ( Z \ 1 „ ( 171 

(!) —(^-~ e 2. 177=) TjM^ 



/ 2n \ 
Vn+m— z/ 

u=0 v v 7 Z=l fc=0 r=l 



4JV+1 / z yZN+l _ -1^, / 2r X ( _ 1)u -fc / m N / m N 2r+fc- U 



/ Tl 



+ 0(e- m2 ' n n- 1 - 2N ) 



as n — > oo. ffere, ffte F r j are some constants the explicit form of which is of no importance in the sequel, and 
Hk(z) denotes the k-th Hermite polynomial, that is, 

m Hk{z) V ( ~ 1)fc ' m (2z)2m ' fc *->o 

W fc! ^ (fe-m)! (2m - jfe)!' " ' 

m>0 v y v y 



The lemma above follows from Stirling's approximation for the factorials. For a detailed proof we refer to [5, 
Lemma 6] . 

3. Height 

In this section we derive asymptotics for the distribution as well as for the moments of the random variable H nyP . 
As mentioned before, the number of p-watermelons with length 2n and height < h is given by = n+v 

Consequently, we have for the distribution of H n ^ p 

(p) 

(3) ¥{Hrwp + l<h} = ^. 
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Theorem 1. For each fixed t € (0, oo) we have the asymptotics 

(4) P <t} = -^i- dot ((-1)^(0) - (t) e-* 2 ) + O (n-^e-r 

as n — > oo, where H a {x) denotes the a-th Hermite polynomial. 

Proof. Set x = (.To, . . . , Xp—i) and y = (yo, ■ ■ ■ , Vp-i)i an d consider the more general quantity 

mW(x,y)= det (Y 2 " ") - f 2 " 

Factoring ( 2 ^) out of each row of the determinant above and replacing each entry with its asymptotic expansions 
as given in Lemma [31 we find the asymptotics 



-L P i(*,y) = (t) P (^(x,y) + O (e- fc, /«„-i-«r)) 



where 

/4N+1 



£> JV (x,y)= det V T tt; ^(0,n)- ^±3 T^.n) 

°- lJ<p V ^ VV V" / \ Vn J 

and iV > is an arbitrary integer. Here, T u] N(h,n) is given by (see LemmaE]) 



/ 17 tU I /~\ 3N+1 U ~ X 21 

T u;N (h, n) = C -*V« M^tM + 5>"« £ J> 

V ' 1=1 fc=0 r=l 



2r \ (-l)"- fc ( h\( h^ 2r+k ~ uS 



hi 



The quantity D^(x, y) is seen to be polynomial in the cc^'s and y/s. This polynomial is divisible by the factors 
(xj — Xi) and (t/j — yi) for < i < j < p, for if Xj = Xi then the j-th and the «-th row are equal and, therefore, the 
determinant is zero (if yj = yi then the j-th and i-th column are equal). Hence, 

n fa - x i)ivj - Vi) 

D N ( x ,y) = n -® °^ <j<P u X (n,h) + 0{n^e^ '">), n - oo. 

n 3- 

0<3<P 

Here, the error term is determined by noting that every power of xj and yj entails a factor of n -1 / 2 , as can be 
seen from the definition of Djv(x, y) above. The unknown coefficient x( n > h) can now be determined by comparing 
coefficients on both sides of the equation above. Comparing the coefficients of IIj=o EjUj' we obtain (after some 
simplifications) the equation 

det f(-iyH i+j (0) - H i+j (4=) e- h2/n ) = h). 



0<i,j<p 

If we specialise by setting Xj = yj = j, then we see that 



■a - »- <B Cr)'^<, (h/«w» - "'+> (£) -'■""■) + (r-'"'-'") ■ 

Setting h = tyfn and replacing rn^l with its asymptotic equivalent as given by Lemma we obtain the result. □ 

Let us now turn our attention to the moments of the distribution of H n p . Clearly, we have for s £ N, 

_„(p) (p) (p) „,(p) 

h>l m ™ h>l m ™ 

The dominant terms of the asymptotics for the moments are going to be expressed by linear combinations of certain 
infinite exponential sums. Asymptotics for these sums are to be determined now. 
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Figure 2. Comparison of the c.d.f. of the random variable "height" on the set of 3— watermelons 
of length 500 without wall (dotted curve) and the limiting distribution as given by Theorem [1] 



Lemma 4. For v > and fi > define 

h>l 

This sum admits the asymptotic series expansion 

W "'~2 I 2 iU + (2to + ^+1)!to! ' 

v ' xrv m>0 v ' 

as n — » oo, where T denotes the gamma function and B m is the m-th Bernoulli number defined via X^>o -^i^Vi' = 
V(e*-1). 

Proof (Sketch). Asymptotics for sums of this form can often be obtain by means of Mellin transform techniques. 
For a detailed overview of Mellin transforms, harmonic sums and asymptotics, we refer to [7]. 
We proceed with a sketch of the proof. The inverse Mellin transform gives 



c+ioo 

dz 



h>l h>\ 
1 

'2ttI 



/,„(») -E^* 

\z) (^f* t{2z-v)d. 



c+ioo 

re 



The integrand has simple poles at z = (v + l)/2 and 2 = 0, — 1, —2, . . . corresponding to the poles of the zeta and 
the gamma function, respectively. The result is now obtained by pushing the line of integration to the left and 
taking into account the residues. 

For the sake of convenience, we mention the evaluations 

Res T(z) = i — m = 0,1,2,..., 

z= — m ml 

ResCOz) = 1 

z — l 

C,{-m)=B m+1 — , m = 0,1,2,..., 

TO + 1 

where B m denotes the m-th Bernoulli number defined via X^>o -^3*Vi' = V ( e * — !)■ ^ 

The rest of this section is devoted to the proof of Theorem [3 below, which gives the final expression for the 
asymptotics of the moments. In order to present the proof of this theorem in a clear fashion we split it into a series 
of lemmas. For a more detailed overview of the proof, we refer directly to the proof of Theorem [21 



6 



THOMAS FEIERL 



As a first step, we prove in Lemma [5] a preliminary asymptotic expression for the moments of the height 
distribution. The presented compact form of the asymptotics makes use of certain linear operators that are going 
to be defined now. 



Definition 1. Let Hi and So denote the linear operators defined by 

Si (h u e-^ h2 ) -L 

Ho = (-D"^T, 

where Bk denotes the A;-th Bernoulli number. 
By Lemma S] we have 



2 V 2 / VM 
B 



- Si (hre-**) n^' 2 + S (Ve"^) + O^" 1 ), 



so that Si and So yield the coefficients of the first two terms in the asymptotic expansion of f v ^{n). 

The preliminary expression for the asymptotics of the moments can now be proven in pretty much the same 
way as in Theorem [T] 

Lemma 5. For s € N, s > 1, the s-th moment of the random variable "height" satisfies the asymptotics 
(6) E (H? l p ) = sSi (k p /i s_1 ) n s l 2 — Hi ( Q K p h s - 2 + r^ 1 ) n^/ 2 + H (k p ) + (n^ 1 ) 

as n — ► oOj where 

Kp = i- det ((-im + ,(o) - (ft) e - ,i2 ) 

1 1 J' 0<t,3<p V / 
0<j<p 

= f ^ 2~(S) /r(-l)^ +j (0)-ff i+j (/ l )e- h2 i/i<p-l\ 

'" '" ' II j!o<j,j< P (-^(Ol-ff^We-* 2 ifi = p-ir 



and 



0<j<p 

Here, H}~(z) denotes the k-th Hermite polynomial. 



Proof. Recall the exact expression for the s-th moment of the random variable "height" (see Equation (|5 

™+2p-2 O) _ (p) 

(7) Efe)= £ (^-(fe-l) s ) " (p) n ' fe - 

h=l m n 

Asymptotics for this quantity can be obtained in pretty much the same way as Theorem [T] Compared to the 
problem of determining asymptotics for the main difference now is the summation over h. 
We consider the more general quantity 

m n ( x , y) - m n,h(x> y) - Q £f <p (( n + Xi - yj) ) ~ o<1 6 /< P ( (n + x 4 - yj) ~ {n + h -Xi- y 3 

where x = (xq, . . . , x p —i) and y = (yo, . . . , y p -i). As a first step, we pull ( ™) out of each row of the determinants 
above. Now, we restrict the range of summation in (0 to 1 < h < n 1 / 2+e for some e > 0. This truncation is 
justified by Stirling's formula, which shows that 



( 2n ) 

\n-j-aJ 



= o 



('--) 



whenever \a\ > n 1 / 2+e . This implies that the total contribution of all summands in ([7]) satisfying h > n 1 / 2+e 
is exponentially small as n — + oo and, therefore, negligible. In all the remaining summands we replace all the 
quotients of binomial coefficients with their asymptotic expansions as given in Lemma [31 Finally, we re-extend 
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the range of summation to h > 1, which, again, introduces an exponentially small error term. This gives the 
asymptotics 

E (Kp) = E [i hS ~( h ~ *)") [^T Dn ^ e) + ° (^ h2/n n^- 1 - 2N )^ , n^oo, 

where e = (0, 1, . . . ,p — 1). Here, the structure of the error term is a consequence of Lemma [5J and the quantity 
Z?at(x, y) is defined by 

(8) zMx,y) = 0< d<* p £ ((«L=£) T u;JV (0,n) 

- J* „ (g' ( ) " w» ») - ) " w, »: 

where A" > is an arbitrary integer and 

K u — k J kl \ \/n 
As a consequence of LemmaSJ we see (after expanding the term (h — l) s ) that 



/ it iu i r~\ 3N+1 u - 1 21 
T U , N (K n) = e-^l- MtLm + ^ n~ l £ £ F r . 

V ' i=l fc=0r=l 



2 ((/I s — {h — If ) O (e-^^n®- 1 -™)) = O (nCD-^+C- 1 )/ 2 ) , 



h>l 

which is negligible for sufficiently large N. Hence, we have the asymptotics 

l2n\V 



( \ 

E (H:j = ±*L. £ ( (h* - (h- 1Y) D N (e, e)) + O (n®"^- 1 )/ 8 ) 

It remains to determine the part of Djv(x, y) that gives the dominant contribution to the asymptotics above. First, 
we note that Djv(x, y) is a polynomial in the Xj's and y,'s. Obviously, -Dtv(x, y) is equal to zero whenever Xj = Xj 
or ?/i = y,j for some i 7^ j, for if Xj = Xj (yi = yj) then the i-th and j-th rows (columns) of the determinants 
involved in the definition of Dat(x, y) are equal, and, therefore, the determinants are equal to zero. This implies 
that Djv(x, y) is of the form 

D N ( X ,y)=n-( ^ i<j<P n % jp ' [x{n,h)+Y^ (^,h)^=+ mM^j + O (tTV* 



0<j<p 



as n — > 00. By comparing coefficients of IIj=o ^jVj on both sides of the equation above, we have already seen (see 
Theorem [l]) that 

-h 2 /n 



X (n, h) = dot ((-1)^(0)) - det ( (-1)^(0) - H l+J ( A 



e 



Analogously we can determine /i). By comparing the coefficients of X& Ilj=o ^j^j on both sides of the equation 
above we obtain the equations 

= £ k (n, h) - Cfe+i(n, h), k < p - 1, 

and 



P"<..i<i» ^(-l)»/f^_,(0) - (^f) e-» '» if i =P- 1, 

Note, that the coefficient of Xk IIj=o x \v\ m ^ n e first determinant of ([5]) is equal to zero, which is easily seen to be 
true for k < p — 1, and for k = p — 1 this is seen to be true by a series of column and row operations that yield a 
new matrix consisting of two non-square blocks. Similar expressions (with i and j interchanged) can be found for 
the rjk(n, h), < k < p. 
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Noting that Hi + j(0) is non-zero if and only if i+j is even we deduce that (— l) l Hi +J (0) = (— l) J f/i +J (0), which 
implies 

£ p -i(n,h) = r) p -i(n,h), 

and also 

p-i 

det ((-im +J (0)) = det (o)) = 2® JJ jl 

3=0 

Here, the last equality has been proven in Lemma [6l 

If we specialise to Xj = yj = j, < j < p, then we obtain 



D 



N 



(e,e)=n"(S) ^(n,ft) + 2^ p _ 1 (n,/ i )n- 1 /^ + O (n^e^ 2 ' n ) , 



where e = (0, 1, ... ,p — 1). 

Choosing N large enough and expanding the term h s — (h— l) s in the asymptotics for E (ff^ ) above, we obtain 
with the help of Lemma 0] the asymptotics 

(2n\P 



and replacing Djv(e,e) with its asymptotic expansion as given above proves the lemma. □ 

Lemma 6. Let Hk(x) denote the k-th Hermite polynomial as defined by Equation We have the determinant 
evaluation 

(9) det ((-1)^/^(0)) =2® jji'- 

Proof. The determinant under consideration is a Hankel determinant. Therefore, we can hope to evaluate it with 
the help of orthogonal polynomials (for details see [15, Section 2.7]). It is well known (see, e.g., [20, page 105]) 
that for k S N we have 

H 2k+1 (0) = and H 2k (Q) = (-l) k ^. 

k\ 

Consequently, we obtain 



det ((-l)^H i+j (0))=2® det r i + ( - 1) ^ 2( ' + 5 2 r^ + ^ + 1 

0<i,j<p \ J 0<i,j<p 



The (i, j)-th entry of the determinant on the right hand side above is seen to be precisely the (i + j)-th moment 
with respect to the Gaussian weight w{x) = -h=eT x ^ on R, that is, 



/2ir 

,-. ; .-x-/2. l + (-l) fc 2 fc /y fk + 1 



x K e~ x /z dx = \ ' =r — — - , k = 0,1, 2 



V27T J-oo 2 y/TT 

The family of monic orthogonal polynomials associated with the weight w(x) is given by 
(10) 2- n / 2 H n (^j, n = 0,1,2,... 

where H n (x) denotes the n-th Hcrmitc polynomial as defined by Equation (J5|). The three term recursion relation 
for the orthogonal polynomials (fT0|) is seen to be (cf. [20, p. 105]) 

2 -,. +1 )/»H„ +1 (£)=x2-*H n (^)-^ (^|). -1,2...., 
with the initial values iJo f^f J = 1 an d 2~ 1 / 2 i?i (^^^) = Now, an application of [15, Theorem 11]) shows that 

o<i,i< P \ 2 0? v 2 ) J -i 

which proves the claim. □ 
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Lemma 7. Let fi > denote a real number. The operator Si from Definition]^ satisfies the relation 
(11) ^(i(h v e-^))={- 1 ifv = ° 

Proof. For v = the claim follows immediately from the definition of the operator Si. For i/>0we calculate 

Si (V+V^ 2 ) = ^Si (h^e-^ 2 ) , 
from which the claims follows upon multiplying by 2/i and rearranging the terms. 

The next result is not obvious at all, and, on the contrary, is a quite surprising fact. 
Lemma 8. Let k p and r p denote the determinants defined in Lemma\S[ We have the relation 



□ 



(12) 

Proof. For the sake of convenience we set 



(p- 1) J^p = t p , p>l. 



fp-i 



J! 



K J=0 

The derivative of a p x p determinant is the sum p determinants, where the j-th addend is equal to the original 
determinant with the j-th row replaced by its derivative. Hence, 



d_ 



>p-2 



where 



/ / #o,o 



Mi = det 



-H l+1 {h)e~ h2 

#i+l,0 



j'=o 

#0,;p-l \ \ 



2— l.p— 1 



-H l+p {h)e- h 



(-l)\ff i+J -(0) -H i+j (h)e 



V V flp-i,o • • • #p-i,p-i / / 

We want to mention that (p— l)CM p _i is equal to the expression for t p except for the constant terms in the last 
row of the determinant. 

For < i < p — 1 the quantity Mi can also be represented by the expression 



Mi = det 



/ / #0,0 



#i-l,0 
#»+l 

(-ly^iWo) 

#i+2,0 



#0,p-l \ \ 



#i-l,p-l 
#i+p 

(-l)' 1+1 ff l+p (0) 

#i+2,p-l 



#p-l,p-l / / 



< i <p- 1, 



V \ # P -i,o 

which is more convenient to work with. 

The Laplace expansion for determinants with respect to the row j + 1, < j < p — 1, gives 

p-i 

Mj = J2(-^y +1 H J+ i+k(0)M^ k , < j < p - 1, 

fc=0 
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where Mj jk denotes the minor of Mj obtained by removing row j + 1 and column k, i.e., 



Mj jfc = dot 



Now, consider the sum 

p-2 



( ( #0,0 



£)i-l,0 



'p-2 p-2 



#0 



fe-1 



fc+1 



£)i-l,fc-l £)i-l,fc+l 
£)i+l,fc-l £)i+l,fc+l 
£)i+2,fc-l £)i+2,fc+l 

£)p-l,fc-l -£)p-l,fc+l 

p-2 



#0,p-l \ \ 



-£)i+l,p-l 
-£)i+2,p-l 



£>p-l,p-l/ / 



E M ^ = EE(- 1 ) j+1 ^+i+fc(°) M ^fc +E(- 1 ) J ' +lfl i+ P (o)M jiP _ 1 . 

The first sum on the right hand side in fact is equal to zero as is going to be shown now. First, note that 

since the matrices involved are transposes of each other. Recalling that H k (0) is non zero if and only if k is an 
even number we deduce that 

(-iy +1 H j+1+k (0)M jtk = -(-l) k+1 H k+1+1 (0)M kj , 

and both expressions correspond to different addends of the double sum above (j + l + k has to be even). This 
shows that the value of the double sum is indeed equal to zero. 
For the second sum we have 

p-2 p-2 
3=0 j=0 

det ( i(- 1 ) kH k+i(0)-H k+ i(h)e- h2 i£k<p-l\ 
o<k% P y\(-l)PH p+l ifk = p-lj' 

which proves the lemma. □ 

We are now able to to state and prove the final expression for the asymptotics of the moments. 

Theorem 2. The expected value of the random variable H n , p satisfies the asymptotics 

(13) E(H n , p )=E 1 (K p )V^ + p-l + 0(n- 1 / 2 y n^^, 
and for s G N, s > 2, we have the asymptotics 

(14) E (H s n p ) = s-!^- V /2 + (« - 1) (p - 1 - |) Si (n p h s - 2 ) n^' 2 + O (n 8 / 2 " 1 ) , n -> oo. 
Here, Hip IS defined by 

k p = 1 - det {{-lfH i+i {Q) H i+j (h) e-» 2 ) , 

0<j<p 

where H k (z) denotes the k-th Hermite polynomial. 

Proof. As a first step we need to establish some simple facts concerning the quantity n p . To be more precise, we 
have to show that n p is an even function with respect to h that has no constant term, i.e., is of the form 

K M 

« P = EE Kmh 2k e- mh2 

k=0 m=\ 

for some numbers K, M and some constants X k m . 
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It is obvious from the definition of the Hermite polynomials (see Equation ) that the k-th Hermite polynomial 
is an even (odd) polynomial whenever k is even (odd). This also implies the equality (— iyHi + j(0) = (— l) J Hi+j(0). 
Now, replacing h by —h in the definition of K p , factoring (— l) 1 out of the i-th row and (— l) 3 out of the j-th row 
we see that the expression remains unaltered. Hence, k p is an even function of h. The constant term of n p is seen 
to be equal to 

'-T^^-det ((-1)^.(0))=!--^-^ ({-l)W'*H i+s (oj)=0, 

where the last equality is a consequence of Lemma El This proves the claimed form of k p . 

We are now going to prove the asymptotics (fl4|) . Therefore, we assume that s > 1. The properties of k p 
established above together with Lemma [7] imply the equation 



and the product rule for the derivative together with Lemma [5] show that 

Si (rph 3 - 1 ) = -(s - - l)Hi ( Kp h s - 2 ) . 

The asymptotics (| 14[) is now obtained from the asymptotics © upon noting that the So-term is negligible for 
s > 2. 

Finally, we prove the asymptotics (|13p and, therefore, assume s = 1. For the sake of simplicity we set 

/p-i N ' 

c=2-& nil 

\3=0 



From Lemma [8] and Lemma [7J we deduce that 

dh'" p J " v ^ V~d/i 



S X (r p ) = (p - l)Si ( = -(p - 1)S X f C-^-xW ) , 



where 

= n det ((-1)% +J (0) - H, l+J (0)e- h2 ) . 

This last determinant can be evaluated to a closed form expression with the help of Lemma El Factoring 1 
(— l) J e~ h out of each column of the determinant we see that 

= ( n (l - (-lj'e"*") ] Q< det ((-lf + ^H i+j (0)) 



yj=0 



hi 



1 - e- 2h ( 1 - e"' 1 



Now, an application of Lemma [7] shows that 



which implies 

s i ( t p) = 1 - V- 

The last step of the proof is the evaluation of the quantity So(k p ). Recalling that k p is an even function with 
respect to h as well as the fact that all odd Bernoulli numbers except for B\ are zero, i.e., B 2v +i = 0, v > 1, we 
deduce the equation 

S ( Kp )-S (l-C X (/ l ))-S h-(l- e - 2 ' 1 ) J. 
The definition of Sq reveals that Sq (h v e~^ h2 ^j is independent of /i. Consequently, we see that 
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Table 1 . This table gives the coefficient of the dominant asymptotic term of KH * as n — > oo 
for small values of s and p (see Theorem [2]) . 



8K^ 


s = 1 


S = 2 


s = 3 


P =l 

p = 2 
p = 3 
p = 4 


±^=0.88... 
a+ 4 ^V5F= 1-51... 
72+45 Vg- 16 ^V^= 1-99. ■■ 

10368+17091^-3776^ nz o on 

20736 — z.oy... 


1 

5 
2 
25 
6 

1915 
324 


fV5r = 1.32... 
20^V5F = 4.45... 
1584+31 3 5 sf" 32v/3 V^ = 9.11... 

520992+165969^-29824^ 1 c A/i 

82944 V n - itKU4 ■ ■ ■ 



This proves the asymptotics (fT3"|) and completes the proof of the theorem. □ 

Table [T] shows the constant of the dominant asymptotic term as n — > oo for the s-th moment of the height 
distribution for small values of s and p. 

4. Range 

We determine the asymptotics for n — > oo of 



1 r 

(15) P{i?„, p <r} = — -y £ Rl+i.,-,,+1 



» 



m « h=2p-2 



Note that m^.j r _/ l+ i — TO i P i r -/i+i ^ s ^ nc number of watermelons with height exactly h and range < r. 
Theorem 3. For each fixed t G (0, oo) we /icwe i/ie asymptotics 

where 

T p (z,w)= o det ((-iy (^H l+J (tz)e-^ 2 ) - (^H l+] {tz + w) e~^ +w A\ . 

- JJ<P V Kiel. J Kiel. J J 

Here, H a denotes the a-th Hermite polynomial. 

Proof. Since rn^\ 2 k = for any k, Equation (|15|) can be rewritten as 



,2p 

P {-Rn,p < r} = n,r ( + ) 1 ' 1 + — { m n:h,r-h+2 ~ m n,h,r-h+l 

The hrst term on the right-hand side is negligible. To see this, we note that m^r+i,! is equal to the number of 
p- watermelons with wall and height < r, which is of order ( 2 ^ l ) P n _p2 as n — * oo (see [5] for details), whereas mffl 

is of order ( 2 ^) P n - ( 2 ) (see Lemma^. 

Asymptotics for the sum on the right-hand side can now be established in a fashion analogous to the proof of 
Theorem [TJ A more detailed presentation of these techniques can also be found in [5, Theorem 2]. We hnd the 
asymptotics 

as n — > oo, where 



T p (t,w) = Q< det ((-IT (j2H i+j (tt)e-^A - 



H l+] {tt + w)e- (lt+w)2 



Now, Taylor series expansion shows that 
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where T" denotes the derivative of T with respect to its first argument. Setting r + 1 = ty/n we see that 



fe=2p-l v ' ^ v v / / fc=2p-l v 



as n — > oo. 



□ 



Remark 2. For the special case p = 1 we recover a well-known fact originally proven by Chung [3] and 
Kennedy [14]. Namely, the equality of the distributions of the height of Brownian excursions and the range 
of Brownian bridges. This result also follows from a more general relation between excursions an bridges proved 
by Vervaat [21]. 

In fact, for p = 1 we have 



-(it+w) 

• - ■_ lei 



which shows that 



2\ e -(ety 



by Theorem [31 This shows that the distribution of the range of 1-watermelons without wall weakly converges to 
the limiting distribution of the height of 1-watermelons with wall (see [5] ) . 
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